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Abstract 



We provide the solutions of the fixed point conditions of the Yukawa sector for a large 
class of N = 1 supersymmetric theories including the minimal and next-to-minimal 
supersymmetric standard models and their grand unified and other extensions. We 
also introduce a test which can discriminate between infra-red stable, infra-red un- 
stable and saddle point solutions, and illustrate our methods with the example of 
the next-to-minimal supersymmetric standard model. We show that in this case, the 
fixed point prediction of the top quark mass is equivalent to that of the minimal 
supersymmetric standard model, supporting previous numerical analyses. 



The fermion mass problem in the standard model arises from the presence of many 
unknown dimensionless Yukawa couplings. A possible solution to this long-standing 
problem is that the standard model could be embedded in some more predictive high 
energy theory with fewer Yukawa couplings. Alternatively one may attempt to relate 
the Yukawa couplings to the gauge coupling, as in the Pendleton- Ross infra-red stable 
fixed point (IRSFP) for the top quark Yukawa coupling |[], or the quasi-fixed point of 
Hill To obtain IRSFPs it is necessary that the unknown dimensionless Yukawa 
couplings are of the same order as the gauge coupling at high energy. Thus at first 
sight this approach would seem to be inapplicable to the small Yukawa couplings of 
the standard model. However string theory yields large Yukawa couplings of order the 
gauge coupling, and this could suggest that the small dimensionless Yukawa couplings 
be reinterpreted as large Yukawa couplings multiplied by a ratio of mass scales. In 
such a scenario fixed points may once again be applicable, as pointed out recently by 
Ross In this case the high energy scale would be the string scale and the infra- 
red region would be near the GUT scale. Thus IRSFPs may provide an unexpected 
resolution of the fermion mass problem. IRSFPs in the MSSM can also be used to 
explain the apparent success of the A T = A& GUT scale relation of some popular 
SUSY GUTs [||. Further fixed points have also been identified in the soft masses 
of supersymmetric theories || [7], §. It is clear that the idea of IRSFPs has general 
applicability, and here we shall follow a model independent approach. 

In this letter we shall present a general approach which may be used to identify 
the fixed points that exist in the dimensionless Yukawa couplings of a large class of 
N = 1 supersymmetric theories, and to determine the stability properties of the fixed 
points in the infra-red. Our assumptions are listed below: 

1. We work at one loop order only. 

2. We assume that the wavefunction anomalous dimensions are diagonal t| oc ^7^. 
This may be due to some symmetry that forbids wavefunction mixing between 
different irreducible representations. For example in the MSSM R-parity forbids 
mixing between the lepton doublet L and the Higgs doublet H\. 

3. We assume a single gauge coupling. This may either be due to several approxi- 
mately equal gauge couplings, which could be the case near the string scale, or 
due to there being a single dominant gauge coupling as in the case of QCD at 
low energies where the electroweak couplings may be neglected. 

Our general approach may be compared to that of Lanzagorta and Ross || for a 
single Yukawa coupling, however unlike these authors we shall present a well defined 
procedure for obtaining the fixed points for theories involving any number of large 
Yukawa couplings in the theory, and having identified the fixed points we shall give 
a general discussion of their infra-red stability. 

1 The relation between the fixed point and quasi- fixed point has been fully elucidated by Lan- 
zagorta and Ross . 
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We shall be interested in the one-loop renormalisation group equations (RGEs) 
for the dimensionless trilinear Yukawa couplings of an N=l SUSY theory, which are 
by now well known and have been widely reported (see for example Martin and 
Vaughn ||). The trilinear part of the superpotential may be written 

W = -h^fafafa, (1) 

where h^ k are the dimensionless Yukawa couplings, fa represent the superfields and 
the sums over i,j, k run over all superfields in the spectrum of the effective theory. 
The relevant RGEs read 

16tt 2 - = h ijp ^ k + ti*> k ji + h pjk jL (2) 

a In// ' F y 

^ 2 4r~ = b a gl (3) 

a in // 



where // is the MS renormalisation scale and the anomalous dimension is defined as 

7*^E h ^ hJm - 25 i9lCt (4) 
z p,i 

to one loop order. b a is the one-loop beta function of the gauge group G a and the 
sum over a here represents a sum over all dominant gauge couplings. C a (i) is the 
quadratic Casimir of the representation of fa under the group with gauge coupling 

9a innt = c a R 8{. (5) 

Here, T$ is a matrix in the R representation for the group labeled by a. Pictorially, 
the RG evolution of the trilinear coupling can be described as an insertion of the 
anomalous dimension correction on each external leg fllQfl . 

According to our second assumption, the wavefunction anomalous dimensions are 
diagonal, 

7i = 1-M (6) 

then the above equations become: 

dh^ k 



where the anomalous dimension is defined as 

1 

p,g 

We now introduce a more convenient notation: 



<h ijk ) 2 n 2 -v- 

16vr 2 ' 16tt 2 ' 16tt 2 ' p v ; 
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in terms of which: 

dY ijk 



dt 



= Y ijk (Ni + Nj + N k ). (10) 



According to our third assumption there is a single gauge coupling in the theory. 
We define the ratio of each Yukawa coupling to the single gauge coupling a as: 

yijk 

R ijk = __ ( U ) 

a 

The RGEs now become: 

d E^ k 

— = a& k (Ni + Nj + N k + ab), (12) 

The non-trivial (non-zero) IRSFPs we seek correspond to dR ^ k = for all k with 
Rijk ^ o, which implies: 

Ni + Nj + N k + ab = 0, (13) 

for all k corresponding to R l ^ k ^ 0. 

In the above notation the anomalous dimension iVj refers to the superfield fa while 
the trilinear coupling between three such superfields is written as a three index tensor 
h l i k . In a given theory each superfield will in general be in some representation of the 
gauge group, and % runs over all gauge indices as well as the different representations 
themselves. Thus most of the entries of the tensor h^ k will be zero due to gauge 
invariance and other imposed symmetries, and many of the non-zero entries will be 
equal by gauge invariance. For example in a toy theory with a single top quark Yukawa 
coupling, assuming that QCD is the only gauge group, the index i would run over 
the 6 components of the quark doublet Q = (t, b), the 3 components of the top field 
t c and the 2 components of the Higgs doublet H, and in general the Yukawa coupling 
would be h l i k which would be a three index tensor with each index taking 11 possible 
values. The top quark Yukawa coupling is characterised by a single parameter h, as is 
clear by dropping all indices and writing the superpotential W = hQt c H 2 . In looking 
for IRSFPs we are interested in the value of h, so it is clear that we must simplify the 
notation in order to remove all the redundant information contained in the tensor h^ k 
to arrive at the physical coupling h. In the case of the top quark Yukawa coupling it 
is clear that we can define Y h = /i 2 /(167r 2 ) and write the RGE for the single physical 
top quark Yukawa coupling as 



dt 

or defining Rh = Yh/a, 

dR h _ 
dt 



dY 

Yh (N Q + N t c + N H2 ) (14) 



R h {N Q + N tc + N H2 +ab). (15) 
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where 



N Q = 2C Q a-Y h 

N t c = 2C t ca-2Y h 

N H2 = -3Y h (16) 

and the QCD Casimirs are Cq = C t c =4/3. The RGE in this case is thus: 

dR h 
dt 

where 



aR h [(r h + b)-QR h ]. (17) 



r h = J2^C R = 16/3. (18) 

R 

and in the MSSM the QCD beta function is b = —3, which leads to the Pendleton- 
Ross fixed point R* h = 7/18. 

In general we can always relabel the Yukawa couplings of any theory as hi where 
i — 1, • • -n now runs over the n non-zero physically distinct couplings in the theory, 
and should not be confused with the previous use of the index i. Defining Yi = 
/t?/(167r 2 ) and R t = Yi/a, we can always write the RG equations for the physically 
distinct Yukawa couplings of the theory as 



dRi ~ r, 

— — = aRi 
dt 



+ 6) - £ S i3 R 3 



(19) 



where rj = J2r 2C_r with the sum being over the representations of the three fields 
associated with the trilinear coupling hi, and SV,- is some matrix whose numerical 
value is fully specified by the wavefunction anomalous dimensions. The fixed point 
is then reached when the right hand side of Eq.|ll] is zero for all i. Thus, writing the 
fixed point solutions as R* there are two fixed point values for each coupling: R* = 
or the solution corresponding to 

[ (r . + 6) - £ S^R*] = (20) 
j 

The non-trivial fixed point solution is: 

n 

^ = E(^%(^ + & ) (21) 

3=1 

We shall begin our discussion by considering only the non-trivial fixed points, and 
later extend our analysis to the general case where some of the couplings are zero at 
the fixed point. 

To determine the infra-red stability of the system in Eq.|T9|, we need to Taylor 



expand it around the fixed point given in Eq.|2l]. We can then drop all except the 



linear terms, the resulting system of which allows an algebraic solution and can thus 
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be tested for infra-red stability. We therefore make a change of variables to p«(t) = 
Ri(t) - R*. The RGE Eq.[l| then becomes 



<]l>i{V] = a(t)(pi(t) + R*) 



(r, + 6 l )-E^(P.( i ) + ^: 



dt 

where we have substituted the fixed point values of R* from Eq.|2T]. 



(22) 



When we drop the quadratic terms in Eq.|22] and change the independent variable 



from t to a by Eq.R we obtain the linearised system 



<IMn ~ (23) 



d\na(t) b 



3=1 



(no sum on i.) Eq|23]then describes the behaviour of the trajectories as they approach 
the fixed point. It has a general solution: 

n 

Pl {t)^Y.^xf ) m)) x \ (24) 

k=l 



.(*) 

equation 



where Xj , X k are the k = 1, ■■■n eigenvectors and eigenvalues of the eigenvalue 



= "^RiSij- (25) 

The 0^ in the solutions are constants which describe the particular linear combination 
of eigenvectors and are set by initial boundary conditions. The infra-red stability 
properties of the solutions in Eq.|23] are independent of the sign of the one-loop gauge 
beta function b, as we shall now show. 

For b > 0, a decreases with decreasing renormalisation scale p. We require every 
eigenvalue of the matrix ^R*Sij to have a positive real part for infra-red stability, 
so that pi — » as p decreases. For b < 0, at increases with decreasing renormalisation 
scale p, in this case we require all the eigenvalues of the matrix — rnR*Sij to have 
negative real parts, which is equivalent to the previous condition for b > 0, so the 
stability condition is independent of the sign of the beta function, as claimed. = 
corresponds to a direction in coupling space which is neither attracted to, nor repelled 
by the fixed point (to first order). For each of these directions there should be 
one free parameter in the solution to the fixed point equations which embodies the 
information of where a solution lies along this direction (and is set by the initial 
boundary conditions). This corresponds to some information about the higher energy 
physics being retained at lower energies. 

Complex eigenvalues always come in complex conjugate pairs, as do their associ- 
ated eigenvectors. Writing X k = + isk, where s k are real, the solution in this 
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case is 



Pi 



E a fc«* j 

k 



l.S k 



(26) 



where the sum runs over all values of k corresponding to pairs of complex eigenvalues. 
Eq.|26| describes a spiral-like solution, the distance to the fixed point being controlled 
by a tk . Thus tk must be positive for the trajectory to be infra-red stable. For b < 0, 
a increases with decreasing normalisation scale so we require Re(Xk) < for all k 
for the solution to be stable in the infra-red direction. If any of the conditions are 
not met, the fixed point is either a saddle point or an ultra-violet fixed point and so 
the fixed point will never be achieved at low energies and must be rejected. In such 
cases we are led to consider other fixed points where some of the couplings take zero 
values at the fixed point. Another motivation for studying such fixed points would 
be if the non-trivial fixed point solution led to values of Ri with unphysical negative 
values (see example later). 

As already remarked, in general for each value of i there are two possible fixed 
point solutions to consider: R* = or its non-zero fixed point value considered 
above. There are thus 2 n fixed points contained in the theory (some of which could 
be degenerate), and so far we have only considered one possibility (the non-zero case 
for each coupling). For each of the remaining possibilities it is straightforward to 
determine the fixed points. Suppose in general that we are considering a possibility 
with m zero fixed point solutions and n — m non-zero solutions. First, we re-order 
the n couplings such that 



[ R i)i=l,...,m 
l^i ] i=m+l,...,n 







E (s-%( rj 

j=m+l 



b) 



(27) 



Note that the non-zero couplings are now determined by the lower right hand (n — 
m) x (n — m) block of the re-ordered matrix S. As we shall see the discussion of the 
non-zero solutions turns out to be similar to that previously considered, while that 
of the zero fixed points is even simpler. 

For the general case including m zero couplings, we need to Taylor expand around 
the fixed point given in Eq.^7]. We make a change of variables to piit) = Riit) — R*, 



where R* = for i 



m. The RGE Eq.fni then becomes 



dpi(t) 



dt 

dpijt) 
dt 



i=l,...,m 



a(t)pi(t) 
a(t)( Pi (t)+R* 



i=m+l,...,n 



&)-E E s tJ ( P ,(t) + r*) 

j=l j=m+l 



+ b)-J2s ijPj (t)- Sij(pj(t) + r*) 

j=l j=m+l 



where we have substituted the fixed point values of R* from Eqf27|. When we drop 
the quadratic terms and change the independent variable from t to a by Eq|| we 



6 



obtain the linearised system 



dpijt) 
dliaa(t) 

dpi(t) 



Pi(t) 



i=l, 



b) 



E Sa R j 

j=m+l 



- i=m+l,. 



1 n 

-R:j2s ljP3 (t) 

° 3=1 



(28) 



dlna(t) 

Eq.|2S| then describes the behaviour of the trajectories as they approach the fixed 
point. 

For the i — 1 



m the equations are of the simple form: 

dpi(t) 



d\na(t) 



\Pi(t) 



- i=l,...,m 



where 



with solution, 



E S^-in + b) 



j=m+l 



\Pi(t)]i= 



l,---m 



Oi(a) 



As in the previous argument, for b > 0, in order to have pi 
require Aj > 0, where here Aj is given by Eq.RO. 



(29) 

(30) 
(31) 

in the infra-red we 



For % = m + 1, . . . , n, assuming that the m zero solutions are infra-red stable, as 
discussed in the above paragraph, we may simplify the procedure by working in the 
infra-red region where [pi]i=i t ... m —> 0, so that we have a simplified set of equations to 
solve: 



dpiit) 
d\na(t) 



1 



(32) 



i=m+l,...,n 



j=m+l 



In fact the argument now parallels that given previously for the non-zero solutions, 
as in Eq.^ 



and the discussion below it, with attention now being focussed on the 
(n — m) x (n — m) lower right hand block of the re-ordered matrix S. 

As a simple example of the application of these results, we consider the example 
of the next-to-minimal supersymmetric standard model (NMSSM) described by the 
superpotential: 

k 

W = hQt c H 2 + XNHi H 2 - -N 3 (33) 



hQt c H 2 + \NH X H 2 - -N 3 

where we have included in addition to the top quark Yukawa coupling the NMSSM 
terms which contain the gauge singlet N. The fixed points of this theory in the 
presence of the gauge couplings have been previously studied numerically JTTJ . It was 
found there that the couplings A and k are attracted towards zero values ||11|| , and it 
is interesting to see how our analytic treatment here reproduces this. The relevant 
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RGEs are JT2], |TI]| , keeping only the QCD gauge coupling: 

dk 



167T" 



16tH 



(iln/x 

din// 
dh 



dlnfi 

In the notation of Eqs.[5], |TT] we find: 

dR k 



k(<ok 2 + 6A 2 ) 
A(2£; 2 + 4A 2 + 3/i 2 ) 
h(X 2 + Qh 2 -2g 2 (C Q + C t c)) 



(34) 



dt 
dR x 

~dT 
dR h 
dt 



aR k [b - QR k - 6R X ] 
aR x [b-2R k -AR x -?>R h \ 
aR h [(r h + b)-R x -6R h } 



(35) 



where r h = 2(C Q + C^). 

We now order the ratios as Ri = (R k , R\, Rh) and write the RGEs for the ratios 
in the form of Eq.|l9| where we identify: 

n = 0, r 2 = 0, r 3 = r h = 16/3 (36) 

with b = —3 for the NMSSM as in the MSSM, and the matrix S given by: 



S 




(37) 



The non-zero couplings are obtained by inverting the matrix S, then using Eq.|21] with 
(rj + b) = (—3, —3, 7/3). The result is: 



R* 



,29 _19 20, 
-18' - 1F' 27' 



(38) 



The negative value of the fixed point for R\ means that we must abandon this fixed 
point physical possibility. 

We next try a fixed point with R\ = 0, and the other two couplings being given 
by their non-zero solutions. Following our procedure we first re-order the couplings 
so that those with zero fixed point values are listed first. In this case we re-order the 
couplings as : Ri = (R\, R k , Rh), with (rj + b) = (—3, —3, 7/3) (unchanged), but now 
we have: 

/ 4 2 3 \ 

S= 6 6 (39) 
V 1 6 / 
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Since by fiat R{ = it only remains to determine the fixed point values for the other 
two couplings. Since the lower 2x2 block of the re-ordered matrix S, is proportional 
to the unit matrix in this case it is easy to invert and we find 

(i?.J?) = (-^) (40) 
For R* k = with the others non-zero, we find 

OK 

<-21-63>- (41) 

The trial fixed point in these cases again contains a non-physical negative value, so 
we abandon these possibilities. 

Finally we try a fixed point with R\ = = 0, but with ^ 0. No further 
re-ordering is required since the two zero ratios occur for i = 1,2. The non-zero fixed 
point for the top Yukawa coupling is given by: 

^ = ^3 = ^33 1 (^3 + &) = ^ (42) 

Having obtained a physically sensible fixed point we must now test its infra-red sta- 
bility. The behaviour of the couplings [Ri]i=i^ around the origin is determined by the 
eigenvalues in Eq.[3(^ which are given by: 

1 25 
[A,.] l=1>2 = -[fiisiJS - + b)] = [-1, --] (43) 

which indicates that this fixed point is attractive in the infra-red direction. The 
behaviour of the coupling R 3 around its fixed point is governed by the single entry 
matrix 

A33 = \r;S 33 = -- (44) 

whose eigenvalue is trivially equal to A = — |. Since b < in this case, a negative 
eigenvalue indicates an IRSFP. In fact R* h = 7/18 is just the usual Pendleton-Ross 
fixed point of the MSSM for the top quark Yukawa coupling discussed previously. 

There are several important technical questions not addressed by the above argu- 



ments. One of these questions is: how fast Jl3[ does the solution approach the fixed 
point? So far, we can only suggest a numerical test where the boundary conditions 
are picked from some range of initial parameters spanning a range centred upon 1 
and of the order 1. We know that the IRSFP is only formally realised in the limit 
fj, — ► 0. While it is clear that when the magnitude of Aj is large, so is the speed of 
the running and therefore the rate at which the fixed point is approached, in models 
of particle physics // never reaches this limit and so the fixed point solution can only 
be used approximately. Some measure on the "error" of this prediction would then 
be useful. It is not clear yet if there are general classes of models which possess an 
IRSFP. Ideally, there would only be one IRSFP in the theory. If two exist in the 
same model, it would be unclear a priori which would be more likely to be realised in 
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nature. Finally, we should note that the solutions R* >m to Eq.^ are not all guaran- 
teed to be positive. Any that are negative will never be reached: the coupling would 
simply approach zero and therefore should be added to the list of couplings i2j=i,..., m . 

To summarise, we have determined the fixed points in a large class of supersym- 
metric models. These fixed points offer the possibility of increased predictivity and 
therefore testability of the models. There are however there are several important 
theoretical questions raised by this approach. Perhaps the most important is that of 
the existence and uniqueness of IRSFPs in any given theory. In the absence of such 
general theorems we have provided a well defined procedure for exploring the fixed 
point properties. Our suggested procedure is the following: start with the case of a 
non-trivial fixed point for all the n physical couplings, and examine its stability; then 
go on to the case with one zero fixed point coupling and n — 1 non-zero fixed point 
couplings and examine its stability; and so on down to the case of the completely 
trivial fixed point with n zero fixed point couplings. We have the given conditions 
for infra-red stability for the general case of m zero fixed point couplings and n — in 
non-zero fixed point couplings, which enables this procedure to be followed. Applying 
the analysis to the NMSSM in the low tan/3 regime, we are able to show that the 
fixed point prediction for the top quark mass is equivalent to the MSSM prediction. 
These techniques may also be applied to the fermion mass problem JTTJ . 
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